r 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-01-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  the  burden  to  Department  of  Defense.  Washington  Headquarters  Services  Directorate  for  Information  Operations  and  Reports 
(0704-0188).  1215  Jefferson  Davis  Highway.  Suite  1204.  Arlington  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law.  no  person  shall  be 
subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

3.  DATES  COVERED  (From  -  To) 

20  Jul  2006 

REPRINT 

4.  TITLE  AND  SUBTITLE 


Gyrotropic  guiding-center  fluid  theory  for  turbulent  inhomogeneous 
magnetized  plasma 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 
61102F 


6.  AUTHORS 

Jasperse,  John  R.,  Basu,  Bamandas,  Lund,  Eric  J.*,  and  Bouhram** 


5d.  PROJECT  NUMBER 
2311 


5e.  TASK  NUMBER 

SD 


5f.  WORK  UNIT  NUMBER 

A3 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  A^SBXP 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL- VS-H  A-TR-2007- 1012 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 
AFRLA^SBXP 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 
Approved  for  Public  Release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 

Reprinted  from  Physics  of  Plasmas, Vo\  13,  No.  7,072903  (2006).  ©  2006,  American  Institute  of  Physics.  DOI:  10.1063/1.2220006 
*  University  of  New  Hampshire,  Durham,  NH  03824  '“'^CETP-CNRS,  4  Avenue  de  Neptune,  94107  Saint-Maur  Cedex,  France 


14.  ABSTRACT 

In  this  paper,  a  new  fluid  theory  is  given  in  the  guiding-center  and  gyrotropic  approximation  which  is  derivable  from  the  Vlasov-Maxwell 
equations.  The  theory  includes  the  effect  of  wave-particle  interactions  for  the  weakly  turbulent,  weakly  inhomogeneous,  nonuniformly  magnetized 
plasma,  and  it  is  applicable  to  a  variety  of  space  and  laboratory  plasmas.  It  is  assumed  that  the  turbulence  in  random  and  electrostatic,  and  that  the 
velocity-space  Fokker-Planck  operator  can  be  used  to  calculate  the  correlation  functions  that  describe  the  wave-particle  interactions.  Conservation 
laws  are  derived  that  relate  the  low-order  velocity  moments  of  the  particle  distributions  to  the  turbulence.  The  theory  is  based  on  the  work  of 
Hubbard  [Proc.  R.  Soc.  London,  Ser.  A  260,  1 14  (1961)]  and  Ichimaru  and  Rosenbluth  [Phys.  Fluids  13,  2778  (1970)].  In  the  work  presented  here, 
the  idea  is  proposed  that  the  fluid  equations  can  be  solved  (1)  by  using  measurements  of  the  low-order  velocity  moments  to  specify  the  initial  and 
boundary  conditions. 


15.  SUBJECT  TERMS 

Vlasov-Maxwell  equations  Fokker-Planck  operator 

Inhomogeneous,  gyrotropic,  magnetized  plasma 


guiding-center 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

C.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

Ramanda.s  Ra.su 

UNCL 

UNCL 

UNCL 

19B.  TELEPHONE  NUMBER  (Include  area  code) 

Standard  Form  298 (Rev.  8/98) 

Prescribed  by  ANSI  Std.  239.18 


AFRL-VS-HA-TR-2007-1 01 2 


PHYSICS  OF  PLASMAS  13,  072903  (2006) 

Gyrotropic  guiding-center  fluid  theory  for  turbulent  inhomogeneous 
magnetized  plasma 

John  R.  Jasperse  and  Bamandas  Basu 

Air  Force  Research  Laboratory,  Space  Vehicles  Directorate,  Hanscom  AFB,  Massachusetts  01731 

Eric  J.  Lund 

Space  Science  Center,  University  of  New  Hampshire,  Durham,  New  Hampshire  03824 

Mehdi  Bouhram 

CEJP-CNRS,  4  Avenue  de  Neptune,  94107  Saint- Maur  Cedex,  France 

(Received  13  March  2006;  accepted  14  June  2006;  published  online  20  July  2006) 

In  this  paper,  a  new  fluid  theory  is  given  in  the  guiding-center  and  gyrotropic  approximation  which 
is  derivable  from  the  Vlasov-Maxwell  equations.  The  theory  includes  the  effect  of  wave-particle 
interactions  for  the  weakly  turbulent,  weakly  inhomogeneous,  nonuniformly  magnetized  plasma, 
and  it  is  applicable  to  a  variety  of  space  and  laboratory  plasmas.  It  is  assumed  that  the  turbulence 
is  random  and  electrostatic,  and  that  the  velocity-space  Fokker-Planck  operator  can  be  used  to 
calculate  the  correlation  functions  that  describe  the  wave-particle  interactions.  Conservation  laws 
are  derived  that  relate  the  low-order  velocity  moments  of  the  particle  distributions  to  the  turbulence. 

The  theory  is  based  on  the  work  of  Hubbard  [Proc.  R.  Soc.  London,  Ser.  A  260,  114  (1961)]  and 
Ichimaru  and  Rosenbluth  [Phys.  Fluids  13,  2778  (1970)].  In  the  work  presented  here,  the  idea  is 
proposed  that  the  fluid  equations  can  be  solved  (1)  by  using  measurements  of  the  turbulence  to 
specify  the  electric-field  fluctuations;  and  (2)  by  using  measurements  of  the  low-order  velocity 
moments  to  specify  the  initial  and  boundary  conditions.  ©  2006  American  Institute  of  Physics. 

[DOI:  10.1063/1.2220006] 


I.  INTRODUCTION 

In  the  literature,  there  have  been  two  traditional  ap¬ 
proaches  to  plasma  turbulence:  one  is  to  view  the  plasma  as 
an  initial-value  problem;  and  the  other,  as  plasma  where  the 
turbulence  is  fully  developed  and  in  a  saturated,  stationary 
state.  In  the  first  approach,  the  plasma  is  initially  in  contact 
with  a  source  of  free  energy  and,  as  a  result,  one  or  more 
collective  modes  of  the  plasma  are  driven  unstable  and  grow 
exponentially  with  time.  The  theoretical  problem  here  is  to 
solve  for  the  time  evolution  of  the  system  of  wave-kinetic 
and  plasma-kinetic  equations.  Quasilinear  theory, more 
general  weak  turbulence  theories,"^’^  and  the  renormalized, 
resonance  broadening  theory^"^  are  successful  examples  of 
this  approach.  However,  if  the  plasma  remains  in  contact 
with  sources  and  sinks  of  free  energy,  the  particle  and  field 
distributions  will  often  be  driven  to  a  turbulent,  steady  or 
quasi  steady  state.  The  plasma  may  be  either  weakly  or 
strongly  turbulent.  In  this  situation,  theories  have  been  de¬ 
veloped  to  describe  the  particle  and  field  distributions  for  the 
fully  developed  turbulent  plasma.  The  main  theoretical  prob¬ 
lem  here  is  to  find  the  correlation  functions  that  determine 
the  anomalous  transport  and  relaxation  processes  that  oper¬ 
ate.  Over  the  years,  treatments  along  these  lines  have  pro¬ 
duced  important  advances^”^^  and  have  been  discussed  and 
extended.*^ 

The  literature  on  anomalous  transport  is  vast.  For  early 
work,  see  Refs.  13-15  and  for  a  review  of  some  applications 
to  the  earth’s  ionosphere,  see  Ref.  16  and  the  citations 
therein.  For  more  recent  work  on  anomalous  transport  in 
axisymmetric,  toroidal  laboratory  plasmas  using  the  drift- 
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kinetic-equation  approach,  see  Ref.  17,  and  for  a  recent  brief 
review  of  applications  to  some  astrophysical  plasmas,  see 
Ref.  18.  Particle  simulations  have  also  been  used  to  study 
anomalous  transport  in  the  earth’s  magnetosphere.’^"^*  For  a 
discussion  of  a  different  approach  to  plasma  turbulence  using 
renormalization-group  methods  and  some  applications  to 
fully  developed  intermittent  turbulence  in  space  plasmas,  see 
Refs.  22  and  23. 

In  this  paper,  we  present  a  new  method  for  incorporating 
plasma  turbulence  into  the  fluid  equations  which  is  appli¬ 
cable  to  a  variety  of  space  and  laboratory  plasmas.  The  fluid 
equations  are  given  in  the  guiding-center  and  gyrotropic  ap¬ 
proximation  for  weakly  inhomogeneous,  nonuniformly  mag¬ 
netized  plasma  where  the  particles  are  transported  in  one 
spatial  dimension  (the  distance  along  the  magnetic  field)  but 
the  turbulence  is  two-dimensional.  In  deriving  them,  we  start 
with  the  Vlasov-Maxwell  equations.""^  We  assume  that  the 
turbulence  is  random  and  electrostatic,  and  that  the  velocity- 
space  Fokker-Planck  operator  can  be  used  to  calculate  the 
correlation  functions  that  describe  the  wave-particle  interac¬ 
tions.  The  method  we  present  is  based  on  the  work  of 
Hubbard,’®  as  well  as  Ichimaru  and  Rosenbluth,”  hereafter 
called  I  and  R.  Since  the  probability  distribution  that  governs 
the  turbulence  may  not  be  Gaussian,  the  Fokker-Planck  op¬ 
erator  is  an  approximation  which  neglects  terms  in  the  ex¬ 
pansion  of  the  distribution  function  in  powers  of  Av  higher 
than  order  2.  The  friction  and  diffusion  coefficients  depend 
on  the  dielectric  screening  function  and  the  spectral  density 
of  the  longitudinal  electric-field  fluctuations  for  the  turbulent 
plasma.  Another  assumption  inherent  in  the  Fokker-Planck 
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method  is  that  the  turbulence  is  sufficiently  weak  so  that  the 
fluctuating  electric  field  has  a  small  effect  on  the  unperturbed 
particle  orbits.  We  also  assume  that  the  dielectric  screening 
function  and  spectral  density  vary  weakly  in  space  and 
slowly  in  time  compared  to  the  strong  spatial  and  fast  tem¬ 
poral  variation  of  the  correlation  function  for  the  electric- 
field  fluctuations.  In  this  way,  the  plasma  may  be  treated  as 
locally  homogeneous  and  stationary.  This  concept  is  similar 
to  that  given  in  Ref.  25:  the  separation  of  space  and  time 
scales  between  (1)  the  turbulence  and  the  fluid  quantities  and 
(2)  the  use  of  a  weak  turbulence  theory  to  treat  the  interac¬ 
tion  between  the  evolving  turbulence  and  the  plasma  par¬ 
ticles.  For  an  interesting  and  immensely  useful  way  of  incor¬ 
porating  kinetic  effects  (Landau  damping)  into  the  fluid 
equations  using  gyrokinetic  theory,  see  Refs.  26  and  27. 

If  we  could  find  a  renormalized  solution  of  the  wave- 
kinetic  and  plasma-kinetic  equations,  the  problem  would  be 
solved.  The  renormalized  propagator  could  then  be  used  to 
find  the  renormalized  dielectric  screening  function  and  the 
renormalized  spectral  density  for  the  longitudinal  electric- 
field  fluctuations.  However,  the  development  of  a  renormal¬ 
ized  solution  for  turbulent,  inhomogeneous,  nonuniformly 
magnetized  plasma  is  indeed  a  formidable  problem.  The  idea 
that  we  present  here  is  to  bypass  this  difficult  problem  by 
solving  the  fluid  equations,  where  measurements  are  used 
to  specify  the  turbulent,  electric-field  fluctuations  and  where 
measurements  of  the  low-order  velocity  moments  of  the 
particle  distributions  are  also  used  to  specify  the  initial  and 
boundary  conditions. 

The  Birkeland  current  system  of  the  earth’s 
magnetosphere^^  is  a  system  of  upward  and  downward  mag¬ 
netic  field-aligned  electrical  currents  that  flow  between  the 
magnetosphere  and  the  ionosphere  at  high  geomagnetic  lati¬ 
tudes.  The  plasma  in  the  earth’s  magnetosphere  is  driven 
unstable  by  its  continuous  interaction  with  the  solar  wind  and 
is  far  from  equilibrium.  The  primary  motivation  for  the  work 
reported  in  this  paper  is  to  develop  a  fluid  theory  applicable 
to  the  Birkeland  current  system  where  the  plasma  is  weakly 

inhomogeneous,  the  geomagnetic  field  is  nonuniform,  and 

28 

electrostatic  plasma  turbulence  is  known  to  occur. 

In  Secs.  II-VI,  we  give  the  fluid  theory  in  the  guiding- 
center  and  gyrotropic  approximation  in  the  presence  of  ran¬ 
dom  electrostatic  turbulence  where  the  particles  are  trans¬ 
ported  in  one  spatial  dimension  (the  distance  along  the 
magnetic  field)  but  the  turbulence  is  two-dimensional.  In 
Secs.  VII  and  VIII,  we  give  the  fluid  equations  for  quiescent 
(nonturbulent),  drifting,  bi-Maxwellian,  electron-ion  plasma 
in  a  nonuniform  B  field  and  show  that  they  may  be  solved  to 
give  the  expected  result  for  equilibrium,  electron-ion  plasma 
in  a  uniform  B  field.  In  Sec.  IX,  we  summarize  and  discuss 
the  results. 


II.  THE  ENSEMBLE-AVERAGED  VLASOV-MAXWELL 
EQUATIONS  FOR  ELECTROSTATIC  TURBULENCE 

The  ensemble-averaged  Vlasov-Maxwell  equations  for 
electrostatic  turbulence  may  be  determined  from  Appendix  A 
by  neglecting  SB  and  choosing  a  time-independent  model  for 
(B):  ^^0,  (B)=B(r).  This  yields 


d 

—  -h  V  • 

dt 


Sr 


(E) 


1 

-h-v  X  B 
c 


(1) 


(2) 


—  •(E)  =  477-p, 

dr 


|;X(E>  =  0. 

dr 


(3) 


Here,  we  can  write  (E) =-<?</>/ (9r,  so  we  obtain  Poisson’s 
equation 


^0=-47rp. 

dr 


(4) 


These  equations  are  given  in  Gaussian  units  and  are  valid  in 
the  standard  six-dimensional  phase  space  discussed  in 
Appendix  A. 

The  conservation  relations  for  C„  for  electrostatic  turbu¬ 
lence  are  obtained  from  (A6)-(A8)  as 

\  dhc„  =  0,  (5) 


S  +  (<5E-)I  =  0, 

„  J  dr  Att  dr  %it 

(6) 

^  j  d^v^myC„+j^:^{SE^)  =  0.  (7) 


III.  THE  KINETIC  EQUATIONS 

IN  THE  GUIDING-CENTER  AND  GYROTROPIC 

APPROXIMATION 

For  problems  where  B  is  a  weakly  varying  function  of  r, 
it  is  reasonable  to  seek  approximate  solutions  of  the  kinetic 
equations  in  the  guiding-center  coordinate  system  which  are 
gyrotropic  at  each  point  along  the  B-field  flux  tube.  For  the 
case  where  the  magnetic  field  is  strong  and  a  weakly  varying 
function  of  position,  and  the  ensemble-averaged  electric  field 
is  weak  and  a  weakly  varying  function  of  position  and  a 
slowly  varying  function  of  time,  the  guiding-center 
approximation^^’^®  can  be  made.  If  Ib±{Ib\\)  denotes  the 
length  scale  for  the  variation  of  B  in  a  direction  perpendicu¬ 
lar  (parallel)  to  B,  then  by  a  weakly  varying  B  we  mean  one 
where  The  lengths  and  are  de¬ 

fined  as 

a„^  =  \vJD.X  a<^i=|ui/n„|,  (8) 

where  Vla=qcfilmcf  is  the  gyrofrequency  of  the  particle  in¬ 
cluding  the  sign  of  q^.  The  length  scales  are  defined  in  the 
standard  way.  For  example,  lB\\=\{B~^dBldsY\  where  s  de¬ 
notes  the  distance  along  B,  with  a  similar  definition  for 
Similar  conditions  for  the  variation  of  (E)  can  be  given.^^  If 
^£i(^£ii)  denotes  the  length  scale  for  the  variation  of  (E)  in  a 
direction  perpendicular  (parallel)  to  (E),  then  by  a  weakly 
varying  (E)  we  mean  one  where  /£i ^ 
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slowly  varying  function  of  time  is  one  that  varies  on  a  time 
scale,  T,  where  r«=|27r/flj. 

The  self-consistent  equations  of  motion  for  the  particles 
in  the  guiding-center  approximation  have  been  found  by 
Sivukhin.^^  Using  Eq.  (12.14)  of  Sivukhin  and  our  notation 
in  (1),  we  have 


1  d 

1  T  +  ''' 

—  +  ^ 

,  V  1 

(E)+  -V  X  B 

1* 

dr  rric^ 

c 

]  (9  d  P.  d  P,  d  ^ 

=  1  -  +  R  •  —  +  — —  +  — —  f  (/■„), 

dv  m^dv\\  rn^dv^] 


(9) 


where  R,  Py,  and  P^  are  given  by  (6.10)-(6.12)  of  Sivukhin. 
If  we  neglect  all  drifts  perpendicular  to  B  (E  X  B,  gradient  B, 
and  curvature  B),  we  note  that  (6.10)-(6.12)  of  Sivukhin  are 


R  =  Uyb, 


(10) 


P,  =  qJE}b-[b-  ,  (11) 


=  [b  •  (<9B/<9r)/2S]m„U|,i;  ^ ,  (12) 


where  b=B/B.  We  now  replace  the/„  dependence  in  by 
the  gyrophase  average  of  /^,  and  take  the  gyrophase  average 
of  (1)  to  obtain 


dt  ”^5  2B  ds 


Sv  I 


A 


=  c 


a» 


(13) 


/a  =  /a('S.'.«^X.l'll)  =  (27r)  '  f  dip{f„), 

Jo 

C„  =  C„(s,Mii,D||)  =  (27r)  '  d(pCa, 

Jo 


(14) 


(15) 


where  s  denotes  the  distance  along  B  and  cp  is  the  gyrophase 
angle,  both  of  which  are  expressed  in  the  guiding-center  co¬ 
ordinate  system.  We  emphasize  here  that  these  equations  are 
valid  only  for  situations  where  transport  along  B  dominates 
transport  perpendicular  to  B.  The  bar  symbol  denotes  the 
gyrophase  average,  which  is  also  calculated  in  the  guiding- 
center  coordinate  system.  The  guiding-center  coordinate 
system^®  is  one  where  the  velocity-space  coordinates  slowly 
change  their  orientation  as  s  varies  so  that  the  axis  is 
always  parallel  or  antiparallel  to  B.  In  general,  the  B  field  is 
curved  and  varies  slowly  in  space  as  described  above.  We 
intQrpTQt  fa(s  J ,v ^ ,V\\)dsv j^dv (^tt)”^  times  the  aver¬ 
age  number  of  particles  per  unit  area  of  type  a  for  which  the 
coordinates  of  the  guiding  center  of  the  motion  lie  between  s 
and  s+ds,  while  the  velocities  and  Uy  lie  between  and 
v^-^dv^  and  un  and  Uh+^/uh,  respectively.  Poisson’s  equation 
becomes 

^(f)lds^  =  -  Airp,  (16) 


where  ((>=(f)(sj)  and 


P  =  p(j,/)  =  2<7„  I  d\f„is,t,V2.,Vn).  (17) 

a  J 

Here,  is  normalized  so  that  id^vfa=n^,  where  is  the 
number  density  and  {E\)=E\\=-d(l)l ds. 

When  the  model  for  B  is  specified  and  when  is  given, 
then  (13)  and  (16)  are  the  system  of  kinetic  equations  in  the 
guiding-center  and  gyrotropic  approximation  that  describes 
the  problem,  subject  to  appropriate  initial  and  boundary  con¬ 
ditions.  When  dBlds  =5^  0,  a  direct  numerical  solution  or  some 
orthonormal  expansion  method  of  solving  ( 1 3)  and  ( 1 6)  can 
be  considered.  For  a  uniform  B  field,  dBlds=0,  and  the  left- 
hand  side  of  (13)  does  not  contain  but  the  right-hand  side 
does.  For  this  case,  we  note  that  an  orthonormal  expansion 
method  on  the  variable  needs  to  be  considered.  However, 
in  the  remainder  of  this  paper,  we  pursue  the  multiconstitu¬ 
ent,  multimoment  fluid  equations  as  a  method  of  solution  for 
the  problem. 

The  conservation  relations  for  in  the  guiding-center 
and  gyrotropic  approximation  are  found  from  (5)-(7)  and  are 

J  d^vC^  =  0,  (18) 

S  f  d^vm^v„C„+y^[(SE:)-(SEj}]  =  0,  (19) 

^  J  dsSn 

2  f  ^/'4/«>l  +  i;^)C„+|A(<5£i)  +  (<5£‘)]  =  0. 

^  J  L  otfSTT 

(20) 

Here,  cylindrical  coordinates  are  implied,  so  jd?v 
=  In  obtaining  these  results,  we  have 

used  the  fact  that  SE^  and  SE\\  are  functions  of  s  and  t  and 
denote  the  perpendicular  and  parallel  parts  of  the  fluctuating 
electric  field,  respectively. 


IV.  THE  MULTICONSTITUENT,  MULTIMOMENT  FLUID 
EQUATIONS  IN  THE  GUIDING-CENTER 
AND  GYROTROPIC  APPROXIMATION 

We  wish  to  calculate  the  multiconstituent,  multimoment 
fluid  equations  in  the  guiding-center  and  gyrotropic  approxi¬ 
mation.  To  do  this,  we  multiply  (13)  by  u"  Uy  and  integrate 
over  velocity  space.^*  We  use  the  following  notation: 

=  (w,/;7„) 

=  2ir  dv^\  dv\iu'*’'v{Jj,s,t,Vx.,v^i), 

J  0  J-00 

(21) 

{v'[,vi-CJ  =  (nMCJ 

=  27rf  dvj^f  ^/I;||l;'L*"l;^|C„(s,/;^;x,l>||), 

Jo  J_oc 

(22) 
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A  =  -B-'dBlds  =  A-'dAlds,  (23) 

where  A  is  the  cross-sectional  area  of  the  flux  tube.  Here,  n 
and  /  are  positive  integers  or  zero.  The  velocity  moments  of 
(13)  are 

-jinJ-Ja)  +  -^(nj  +  1  ;7J  +  A(1  +n/2){n,l  +  1  ;/„) 
at  as 

+  l{(qjmjd(f>/ds(nj-  1  ;7j  -  {AH) 

X(M  +  2,/-l;7j}  =  (n,/;Cj.  (24) 

The  lowest  four  velocity  moments  [(«,/) 

=  (0,0),(0, 1),(0,2),(2,0)]  and  Poisson’s  equation  are 

+  B^{n^JB)  =  0,  (25) 

dt  3s 

d  d  IdB  (  \  d(f) 

-  i  *  "“r-  ■  * 

=  Wdi.  (26) 

d  d  1  dB  d(f) 

flaQo^\~  .  ^a{QoA\~  ^al)  Q(/^a^a  ^ 

dt  ds  B  ds  ds 

=  VVdi,  (27) 

,  (28) 

dt  ds 

^  =  -  4ir2  (29) 

as  p 

Here,  we  have  defined  the  velocity  moments  using  the 
following  notation:  n„=(0,0;7j,  (0, 1  ;7a), 

=  (mJ2)(0,2-Jj,  /i„H'„x  =  (m„/2)(2,0;7„),  «a9aii=("'a/2) 
X(0,3;7<»).  and  nWai  =  ('««/2)(2, 1  ;7„).  For  particles  of 
type  a,  is  the  number  density,  is  the  parallel  drift 
velocity,  w^|  is  the  total  parallel  energy  per  particle,  is 

the  total  perpendicular  energy  per  particle,  is  the  total 
parallel  energy  flux  per  particle,  and  is  the  total  perpen¬ 
dicular  energy  flux  per  particle.  By  total,  we  mean  the  sum  of 
the  drift  and  random  parts.  Note  that  is  the 

total  parallel  (perpendicular)  energy  flux  and  /iaW^|(n^w^^) 
is  the  total  parallel  (perpendicular)  energy  density.  The  defi¬ 
nitions  for  q^ii  and  q^x  should  not  be  confused  with  the  heat 
flux  per  particle,  also  denoted  by  q^^  and  <7^  in  Braginskii,^^ 
but  defined  differently  here.  Also,  we  have  changed  the  no¬ 
tation  used  for  the  energy  fluxes  per  particle  from  what  was 
used  in  Ref.  3 1 .  In  our  new  notation,  the  lowercase  symbols 
are  used  for  per  particle  quantities  and  uppercase  symbols 
are  reserved  for  per  unit  volume  quantities. 

We  have  also  introduced  the  following  notation: 

=  m,(0,ljCj,  W^,=(m,/2)(0,2;CJ,  and  W«^  =  (m,/2) 
X(2,0;CJ.  Here,  M^w  is  the  rate  of  transfer  of  momentum 
per  unit  volume  for  particles  of  type  a  due  to  wave-particle 
interactions,  and  Waii(W^i)  is  the  rate  of  transfer  of  parallel 


(perpendicular)  energy  per  unit  volume  for  particles  of  type 
a  due  to  wave-particle  interactions.  The  transfer  rates  are 
functions  of  s  and  t.  The  are  related  to  the  anomalous 
(turbulent)  resistivity  for  the  problem,  and  and  W^x  ^^e 
the  anomalous  (turbulent)  parallel  and  perpendicular  heating 
or  cooling  rates  per  unit  volume  for  particles  of  type  a.  If  the 
momentum  (or  energy)  transfer  rate  is  positive,  then  momen¬ 
tum  (or  energy)  is  gained  by  the  particles.  If  it  is  negative, 
then  momentum  (or  energy)  is  lost  by  the  particles.  These 
quantities  are  not  independent  but  are  related  to  the  turbulent 
fluctuations  through  the  conservation  relations  given  by  (19) 
and  (20). 

Equations  (25)-(29)  are  a  set  of  4N+ 1  equations  for  the 
6A^+1  unknowns:  where  N 

is  the  number  of  plasma  constituents.  Instead  of  using  this 
set,  we  may  introduce  the  parallel  and  perpendicular  tem¬ 
peratures  in  energy  units^^  defined  as 


nJ'J2  =  pJ2  =  {mJ2)  J  d^v(v„  -  m  J Va. 

(30) 

i=Pai  =  {mJ2)  J  d^vv\f„. 

(31) 

It  follows  that  Wa\\=Ta\\/2+maU^J2  and  Wax  =  T'ai 
ing  for  Wa\\  and  in  (26)-(28)  gives 

.  Substitut- 

d  d 

—n„  +  B—{n„uJB)  =  Q, 
dt  ds 

d  d  ^ 

-/njiaUa  +  J  -  d  TT'*" 

dt  ds  B  as 

(32) 

2  d(f) 

^  (7'all  "I"  n  ~ 

ds 

(33) 

\  d  ^  2,  d  \  dB  ^ 

--'iJT'aii  +  +  —riaqaw "  “T" “ 

Z  dt  ds  B  ds 

dcj) 

-H  qoZ^cMa  ^  ■” 
ds 

(34) 

d  d 

dt  ds 

(35) 

1 

II 

(36) 

Here  (32)-(36)  are  a  set  of  4A^+ 1  equations  for  the  1 

unknowns:  7^,,  and  <p. 

The  fluid  equations  for  turbulent,  inhomogeneous,  non- 
uniformly  magnetized  plasma  in  the  guiding-center  and  gy- 
rotropic  approximation  are  given  in  two  representations  by 
(25)-(29)  and  by  (32)-(36).  They  are  the  basis  of  our  fluid 
theory.  For  both  sets  of  equations,  2N  closure  conditions  are 
needed.  When  the  model  for  B  is  given  and  when  Ma\\y^a\\ 
and  Wax  specified,  then  either  set  of  fluid  equations  may 
be  solved  subject  to  appropriate  initial  and  boundary  condi¬ 
tions  to  give  the  weak  spatial  and  slow  temporal  evolution  of 
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the  turbulent  system.  Examples  of  how  this  is  done  for  qui¬ 
escent  plasma  are  given  in  Secs.  VII  and  VIII. 

It  is  interesting  to  compare  the  above  equations  to  the 
standard  approach  for  a  uniform  B  field.  We  see  that  (32)  is 
the  one-dimensional  equivalent  of  (1.11)  of  Braginskii.^^  We 
also  see  that  (33)  is  the  one-dimensional  equivalent  of  the 
left-hand  side  of  Eq.  (1.12)  of  Braginskii,  where  no  drifts 
perpendicular  to  B  or  off-diagonal  elements  of  the  pressure 
tensor  appear.  When  we  add  (34)  and  (35)  for  a  uniform  B 
field,  we  see  that  the  sum  is  also  the  one-dimensional  equiva¬ 
lent  of  the  left-hand  side  of  Eq.  (1.13)  of  Braginskii. 


V.  CONSERVATION  LAWS  IN  THE  GUIDING-CENTER 
AND  GYROTROPIC  APPROXIMATION 

The  conservation  laws  for  the  number  density,  charge 
density,  momentum  density,  and  total  energy  density  for 
electrostatic  turbulence  in  the  guiding-center  and  gyrotropic 
approximation  may  be  written  from  the  results  in  Secs.  Ill 
and  IV.  The  conservation  of  number  density  is 

d  d 

—n„  +  B—{njiJB)  =  Q.  (37) 

dt  Ss 

The  conservation  law  for  the  charge  density  follows  from 
(37)  and  is 

^p  +  B^{j,/B)  =  0,  (38) 

at  ds 

where  we  have  defined  From  (19)  and  (33),  we 

a 

obtain  the  conservation  law  for  momentum  density 

d  d  2\  ^ 

-2j  +  -r2j  7 TT 

at  ^  as  ^  B  as 


X  2  +  mjul  -  -  (<5£?)] 

„  OStSTT 

=  p£„.  (39) 

In  order  to  obtain  the  conservation  law  for  the  total  energy 
density,  we  add  (34)  and  (35)  and  use  (20)  to  obtain 

“  ~~T  2  +  ^al)  + 

B  as  ^  atoTT 


(40) 

For  B  equal  to  a  constant,  (37),  (39),  and  (40)  reduce  to  the 
one-dimensional  forms  of  ( A9)-(A  1 1 )  for  electrostatic  turbu¬ 
lence,  given  in  Appendix  A.  Also  note  that  (37),  (39),  and 
(40)  reduce  to  the  form  given  by  Davidson,^  Chap.  8,  if  B  is 
uniform,  the  plasma  is  homogeneous  and  one-dimensional, 
and  the  turbulence  obeys  the  quasilinear  assumptions. 

It  is  interesting  to  examine  the  integrated,  total  energy- 
density  conservation  equation  for  steady-state  conditions.  In¬ 
tegrating  (40)  from  s^  to  ^2,  we  obtain 


r^2  ] 

+  ds-—'2.nMM  +  (laJ 

a  '  Js,  Bds  „ 


(41) 


Let  US  assume  that  both  q^w  and  are  positive  (outgoing 
particles)  on  the  flux  tube  from  Si  to  S2.  If  B~^dB/ds  is  nega¬ 
tive,  and  if  is  positive  on  the  flux  tube,  y||£||  provides  the 
rate  of  change  of  energy  density  to  drive  the  two  terms  on  the 
left-hand  side  of  (41).  Further,  for  the  case  of  a  uniform  B 
field,  then 


a 


dsjnEfi. 


(42) 


This  means  that  the  total  energy  flux  exiting  the  flux  tube  at 
^2  is  greater  than  that  entering  at  ^i,  if  7||£||  is  positive  on  the 
flux  tube. 

The  conservation  laws  for  turbulent,  inhomogeneous, 
nonuniformly  magnetized  plasma  in  the  guiding-center  and 
gyrotropic  approximation  are  given  by  (37)-(40).  Since  all 
drifts  perpendicular  to  B  have  been  neglected,  the  particles 
are  responding  in  one  spatial  dimension,  s,  but  the  electro¬ 
static  turbulence  is  two-dimensional,  as  (^^}=(^£'^) 
+{SE^). 


VI.  CALCULATION  OF  THE  MOMENTUM 
AND  ENERGY  TRANSFER  RATES 
IN  THE  GUIDING-CENTER  AND  GYROTROPIC 
APPROXIMATION 

The  wave-particle  transfer  rates  per  unit  volume  were 
defined  in  Sec.  IV  as 


0 

1 

VV,,, 

rriaplll 

mav\l2 

We  now  assume  that  the  fluctuating  electric  field  is  random 
(Markovian),  that  the  length  and  time  scales  for  the  one- 
particle  distribution  functions  and  the  two-particle  correla¬ 
tion  functions  separate,  and  that  is  given  by  a  velocity- 
space  Fokker-Planck  operator,^^’^^ 

C„  =  -^-(F{  +  F^)/„+~|;:iyj„,  (44) 

ay  2  as  ay 

where  F^,  F^,  and  are  functionals  of  e  and  {\SE^\),  and 
where  e  and  are  the  dielectric  screening  function  and 

spectral  density  of  the  longitudinal  electric-field  fluctuations 
for  the  turbulent  plasma,  respectively.  By  separation  of 
length  and  time  scales,  we  mean  that  there  is  a  length  /  and 
time  T  such  that  l\>l>l2^^D€  ^  ^2»  where 

the  electron  Debye  length  and  T2  is  discussed  in  Appendix  B. 
Here,  /j  and  Tj  are  the  characteristic  length  and  time  scales 
for  the  one-particle  distribution  functions  and  their  moments, 
respectively.  It  is  important  to  point  out  here  that  these  scal¬ 
ing  assumptions  render  (44)  for  approximate:  the  Fokker- 
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Planck  method  is  valid  only  for  large  length-  and  time-scale 
separations. 

In  the  guiding-center  and  gyrotropic  approximation,  we 
consider  a  Cartesian,  velocity-space  coordinate  system  which 
slowly  changes  its  orientation  as  s  varies  along  B  so  that  the 
axis  is  always  parallel  or  anti  parallel  to  B.  See  Sec.  HI. 
Thus,  vl=v\  and  v\=v\-\-v^y.  Using  integration  by  parts  and 
the  asymptotic  properties  of  for  large  v,  we  see  from  (43) 
and  (44)  that 

=  d^v(27r)~'j  d<p(F{+ (45) 


Waii  =  '^aj  'j  d(p{v^(F{+F^)  +  (l/2)DQf„, 


j  d^v(2-n-)  '  j  d(f>{vj,(F{+ F^) 

+  Vy(F{  +  F?)  +  ( 1/2)(Z>^  +  Z>J  J}/„. 


(47) 


(46) 


In  (45)-(47),  there  is  an  explicit  dependence  on  /„  as  shown, 
and  an  implicit  dependence  on  in  F^,  F^,  and  through 
the  Vlasov-Maxwell  hierarchy  of  kinetic  equations.  As  dis¬ 
cussed  in  Sec.  Ill,  we  note  that  the  guiding-center  and  gyro- 
tropic  approximation  is  the  one  obtained  by  replacing  both 
dependences  by  /^.  In  Appendix  B,  we  calculate  F{,  F^,  and 
in  the  guiding-center  and  gyrotropic  approximation.  We 
obtain  and  as  generalizations  of  Eqs.  (64)  and  (63) 
of  Ichimaru  and  Rosenbluth.^^  Using  the  I  and  R  formalism, 
we  may  also  obtain  an  expression  for  From  Appendix 
B,  we  obtain 


•)> 


(48) 


X  j^(|^E^|(i,/;k,  w))  +  /:||i;||^(|^E^|(5,r;k,  w))  j  +  j  J^(^„)4»j„  Im[e(i,/;k,  w)‘']|  <5(Mfl„  +  AjuUu  -  w),  (49) 


+  ^pj/ifl^‘(^J^(|^£^|(5,t:k,w))  +  ^^jnfl^^(^„)4ffi„Ini[e(s,r;k,w)’']  |^nn„  +  A:||t;||- w),  (50) 


where  to  include  the  sign 

of  Also,  y„(^)  is  the  usual  Bessel  function  of  order  n, 
kj^(kii)  is  the  perpendicular  (parallel)  part  of  k,  and  k}=k\ 
+kl.  The  other  symbols  have  their  usual  meaning.  As  dis¬ 
cussed  in  Appendix  B,  the  transfer  rates  are  weak  functions 
of  5  and  slow  functions  of  r,  as  are  the  low-order  velocity 
moments  of  the  one-particle  distribution  functions.  Here,  we 
reiterate  from  Appendix  B  that  the  following  symmetries 
hold:  (|^£^|(5,r;k,a>)}=(|<5£^|(5,r;-k,-a)));  and 

8(^,/;k,a))*  =  e(5,r;-k,-a)). 

In  deriving  (48)-(50),  we  have  neglected  spatial  diffu¬ 


sion  across  B  and  assumed  that  velocity-space  transport 
along  B  dominates  spatial  transport  along  B.  In  order  to  es¬ 
timate  the  magnitude  of  the  spatial  diffusion  processes  to  see 
if  they  are  negligible  for  a  given  problem,  the  formulas  given 
by  Ichimaru  and  Tange^"^  and  Tange^^  may  be  used.  We  also 
note  that,  since  B  is  a  weak  function  of  5,  so  are 

The  momentum  and  energy  transfer  rates  per  unit  vol¬ 
ume  for  turbulent,  inhomogeneous,  nonuniformly  magne¬ 
tized  plasma  in  the  guiding-center  and  gyrotropic  approxima¬ 
tion  are  given  by  (48)-(50).  They  are  generalizations  of  the 
formulas  given  by  I  and  R. 
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VII.  THE  FLUID  EQUATIONS  FOR  QUIESCENT, 
DRIFTING,  BI-MAXWELLIAN, 

ELECTRON-ION  PLASMA  IN  A  NONUNIFORM 
B  FIELD 

In  order  to  find  the  fluid  equations  for  quiescent,  drifting, 
bi-Maxwellian,  electron-ion  plasma  in  a  nonuniform  B  field 
in  the  guiding-center  and  gyrotropic  approximation,  we  must 
evaluate  the  momentum  and  energy  transfer  rates  due  to 
wave-particle  interactions.  Let  us  assume  that  is  a  drifting 
bi-Maxwellian,  and  that  the  relative  drift  between  the  elec¬ 
trons  and  ions  is  sufficiently  small  so  that  the  plasma  remains 
stable  (quiescent)  as  it  evolves  in  space  and  time.  The  drift¬ 
ing  bi-Maxwellian  distributions  are 


1 


Xexp 


(51) 


where  Va\\,  and  Va±  depend  on  s  and  f,  and  we  have 

introduced  the  thermal  velocities  for  the  particle  distributions 
which  are  given  by  and  vl^  =  Ta±/ma.  We  note 

that  fa  has  been  normalized  to  Substituting  (51)  into 
(48)-(50),  we  see  that  the  velocity-space  integrals  are 
known.^^  We  obtain 


+  ( w  -  nil„  -  kifUj 


+  47^1  Im[e(i,r;k,a))q^]  f  A„(/3„)exp|  - 


{(o-nClg-ktiuJ 

2ky„„ 


(52) 


y-nCl^ 


+  ((<>-  nH„  -  k^uj 


(|^£'‘|(5,/;k,  <!>)>, 


qu 


+  47,^1  Ini[e(j,r;k,w)^]  )A„(/3„)exp|  - 

'  '  2*1, 


(53) 


+  (w  -  nO.„  -  kfiuj 


<|^£^|(5,;;k,w)> 


qu 


+  47„ii  Im[e(s,f;k,w)q^]  [A„(;8jexp|  - 


{(o-nCl„-knuy 


(54) 


r 


Here,  we  have  defined  A„(jc)=exp(-x)/„(x),  ^a=k\v\^/H\, 
and  Inix)  is  the  Bessel  function  of  imaginary  argument  of 
order  n.  The  subscript  “qu”  refers  to  quiescent  plasma.  In 
obtaining  (52)-(54),  we  have  used  some  well-known 
integrals^^  involving  j\{^)  and  the  recursion  properties'^  of 

/.(a. 

In  what  follows,  we  give  explicit  expressions  for  Eqy 
and  (I^E^Dqu  as  weak  functions  of  s  and  slow  functions  of  t 
by  making  use  of  the  separation  of  spatial  and  temporal 
scales  between  the  one-particle  distribution  functions  and  the 
two-particle  correlation  functions.  The  expression  for  the 


quiescent  dielectric  screening  function  for  this  case  is  given 
by 

c(5,  t ,  k,  0>)qy 


=  1+  2  ^  1+  2 


<  I  ^ 


1  + 


7aii  j  riQ.„ 

T„,  /  (I)- nCl„-k«u„ 


W\ 


>-nn„-k„u. 


- 1 


A.(y8„)  , 


(55) 


where  Mz)  is  the  plasma  dispersion  function'^  given  by 
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W(2)  =  \  -z  exp(-  z^/2)  j  dyexpiy^l)  -h  i(7r/2)^^^z 

Jo 

Xexp(-z^/2),  (56) 

and  where  we  have  defined  The  spectral 

density  of  the  longitudinal  electric-field  fluctuations  can  be 
calculated  using  the  superposition  technique  for  the  dressed 
test  particles. For  a  drifting  bi-Max wellian  with 
=  g,  the  spectral  density  is 


(|5£:-|(s,r;k,a))>q„  = 


5^(s,r;k,aj)qu  = 


(57) 


2  AM 


tions  where  electron-  and  ion-neutral  particle  collisions  are 
negligible  compared  to  charged-particle  collisions.  Such  an 
application  is  currently  under  study. 


VIII.  THE  SOLUTION  OF  THE  FLUID  EQUATIONS 
FOR  EQUILIBRIUM  ELECTRON-ION  PLASMA 
IN  A  UNIFORM  B  FIELD 

In  this  section,  we  give  the  solution  of  the  fluid  equa¬ 
tions  for  equilibrium,  electron-ion  plasma  in  a  uniform  B 
field  and  obtain  the  expected  result.  We  consider  the  special 
case  when  the  plasma  is  in  thermal  equilibrium  and  the 
B  field  is  uniform,  i.e.,  when  /z^=n/=Ai,M^=w,=0,T’^|i=7'^^ 
=  r,|i=7’,^  =  7’.  From  (51),  we  see  that 


1 


Xexp 


(b>-nO,„-k^^uy- 


2ttv 


qu 

\3/2 

2  exp 


2vl 


(61) 


(58) 


Since  and  (|<S£'^|)qu  are  known  functions  of  /i^,  m„, 

and  Tal,  then  and  are  also  known  functions 

of  T„\\,  and  7^^,  which  are,  in  turn,  functions  of  s  and 

t. 

As  is  the  case  with  all  moment  equation  descriptions, 
closure  conditions  must  be  specified.  For  electron-ion  plasma 
with  drifting,  bi-Maxwellian  distributions,  the  2A(=4)  clo¬ 
sure  conditions  are  found  by  taking  the  appropriate  velocity 


and  that  (52)-(58)  pass  to  their  equilibrium  forms.  Here, 
v\=Tlma.  In  particular,  the  equilibrium  forms  of  and 
(l^f^Dqu  are 


e(k,  a>)gg  =1+2  . 2  I  1  +  2 

ce=e,i  ^  I  n=-« 


(o  -  nCl. 


K(/Ja) 


X 


w 


)-nCir 


-  1 


(62) 


moments  of  (51).  The  expressions  for  the  perpendicular  and 
parallel  total  energy  fluxes  per  particle  in  terms  of  the  low- 

(|5£2|(k,a,)),,- 

order  velocity  moments  are 

(59) 

j  /!=+» 

^  2  \1/2|I,  1  ^  A„(^a) 

a^e.i  (2'n-)‘'1<:|||f„„=.oc 

^aii  =  “a{(3/2)7^|  +  (l/2)m^l}. 

(60) 

(w-nOJ^ 

-  -,,,2.2 

(63) 

The  4A+1(=9)  fluid  equations  given  by  (32)-(36)  and 
the  2A(=4)  closure  equations  given  by  (59),  (60),  and  (52)- 
(54)  for  the  momentum  and  energy  transfer  rates  when 
supplemented  by  (55)-(58)  represent  a  closed  set  of  self- 
consistent  equations  for  quiescent,  drifting,  bi-Maxwellian, 
electron-ion  plasma  in  a  nonuniform  B  field  in  the  guiding- 
center  and  gyrotropic  approximation.  There  are  6A-i-l(=13) 
equations  for  the  6A-i-l(=13)  low-order  velocity  moments 
and  <f>.  The  equations  may  be  solved  subject  to  appropriate 
initial  and  boundary  conditions  to  give  the  weak  spatial  and 
slow  temporal  evolution  of  the  system. 

We  also  note  that,  since  the  magnetized  Balescu-Lenard 
collision  operators  for  the  multiconstituent  plasma  can  be 
recast  into  Fokker-Planck  forms  (see  Ref.  12),  the  effect  of 
charged-particle  collisions  is  contained  in  the  expressions  for 
the  transfer  rates  given  by  (52)-(58)  for  quiescent,  magne¬ 
tized,  bi-Maxwellian  plasma.  For  an  explicit  solution  of  the 
dielectric  function  for  the  unmagnetized,  Balescu-Lenard- 
Poisson  equations,  see  Refs.  39  and  40.  When  gravity  is 
included  in  the  theory,  (32)-(36)  and  (52)-(60)  may  be  ap¬ 
plied  to  the  earth’s  topside  ionosphere  for  quiescent  condi- 


Using  (56)  for  the  definition  of  the  W  function,  we  find  from 
(62)  and  (63)  that 


Im 


0)  /  Le(k,C0)eqJ’ 


(64) 


which  is  the  relationship  that  is  required  by  the  fluctuation- 
dissipation  theorem."^* As  a  consequence  of  (64),  we  find 
from  (52)-(54)  that  W^^=0.  From  (32)-(36),  we 

see  that  the  equilibrium  solution  for  the  fluid  equations  is 


rig  =  ni  =  n  (a  constant) , 

T'eW  =  Tg^  =  7,11  =  =  7  (a  constant), 

Ug  =  M,  =  0, 

£ii  =  0, 


(65) 

(66) 

(67) 

(68) 
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««ii  =  9<'i  =  ^'iil  =  9a  =  0.  (69) 

The  solution  is  the  expected  result.  It  is  important  to  note 
that  the  fluid  equations  obey  the  fluctuation-dissipation  theo¬ 
rem  in  the  limit  as  the  plasma  approaches  thermal  equilib¬ 
rium.  We  also  note  that  this  result  is  independent  of  the  clo¬ 
sure  assumptions  as  qe\\=(iei-^i\\-QiL-^  symmetry,  if 

Ug=Ui=0. 

IX.  SUMMARY  AND  DISCUSSION 

This  work  was  motivated  by  the  success  of  the  I  and  R 
method  in  treating  temperature  relaxation  problems  for  ho¬ 
mogeneous  plasmas  in  a  uniform  magnetic  field. We  have 
generalized  their  work  in  order  to  obtain  a  fluid  theory  ap¬ 
plicable  to  space  and  laboratory  plasmas,  when  random  elec¬ 
trostatic  turbulence  is  present,  the  magnetic  field  is  nonuni¬ 
form,  the  spatial  gradients  are  weak,  the  guiding-center  and 
gyrotropic  approximation  is  generally  valid,  and  transport 
parallel  to  B  dominates  transport  perpendicular  to  B. 

The  main  results  of  this  paper  are  summarized  as  fol¬ 
lows: 

(1)  We  have  derived  a  set  of  kinetic  equations  in  the 
guiding-center  and  gyrotropic  approximation  where  all 
the  drifts  perpendicular  to  B  are  neglected.  These  equa¬ 
tions  are  valid  for  a  weakly  inhomogeneous  magnetic 
field,  provided  that  the  conditions  on  B  and  (E)  dis¬ 
cussed  in  Sec.  Ill  are  satisfied,  that  the  turbulence  is 
electrostatic,  and  that  the  (E^)xB  drift  is  small  com¬ 
pared  to  the  electron  and  ion  drifts  along  B.  See 
(13)-(15)  and  the  self-consistent  Poisson’s  equation 
given  by  (16).  The  kinetic  equations  include  the  effect  of 
wave-particle  interactions  due  to  the  turbulence. 

(2)  Two  sets  of  multiconstituent,  multimoment  fluid  equa¬ 
tions  valid  in  the  guiding-center  and  gyrotropic  approxi¬ 
mation  are  given  by  (25)-(29)  and  (32)-(36).  The  fluid 
equations  are  new  and  include  anomalous  transport  of 
momentum  and  energy  due  to  the  electrostatic  turbu¬ 
lence. 

(3)  Conservation  laws  are  also  given  in  the  guiding-center 
and  gyrotropic  approximation  that  relate  the  spectral 
density  of  the  longitudinal  electric-field  fluctuations  to 
the  low-order  velocity  moments  of  the  one-particle  dis¬ 
tribution  functions.  See  (37)-(40).  The  conservation 
laws  are  also  new. 

(4)  If  we  assume  that  the  turbulent  electric  field  is  random 
(Markovian),  that  the  length  and  time  scales  between  the 
one-particle  distribution  functions  and  the  two-particle 
correlation  functions  separate,  and  that  the  correlation 
functions  are  given  by  the  velocity-space  Fokker-Planck 
operator,  then  explicit  expressions  for  the  wave-particle 
transfer  rates  in  the  guiding-center  and  gyrotropic  ap¬ 
proximation  can  be  found.  We  have  assumed  that  the 
turbulence  is  sufficiently  weak  so  that  the  fluctuating 
electric  field  has  a  small  effect  on  the  unperturbed  par¬ 
ticle  orbits.  The  transfer  rates  are  functionals  of  e  and 
(|<5£^|},  the  dielectric  screening  function,  and  the  spec¬ 
tral  density  of  the  longitudinal  electric-field  fluctuations 
for  the  turbulent  plasma,  respectively.  Here,  we  have 


generalized  the  work  of  I  and  R  by  introducing  the  ten- 
sorial  spectral  density  of  the  longitudinal  electric-field 
fluctuations  in  the  guiding-center  and  gyrotropic  ap¬ 
proximation.  See  (BIO).  The  transfer  rates  are  given  by 
(48)-(50).  If  large-amplitude  coherent  structures  are 
present  in  the  plasma  so  that  the  length  and  time  scales 
do  not  separate,  then  the  method  presented  here  may  not 
apply. 

(5)  For  quiescent,  drifting,  bi-Maxwellian,  electron-ion 
plasma  in  a  nonuniform  B  field,  we  give  explicit  expres¬ 
sions  for  the  momentum  and  energy  transfer  rates.  See 
(52)-(58).  In  the  limit  as  the  plasma  approaches  thermal 
equilibrium,  the  Fokker-Planck  expressions  for  the  cor¬ 
relation  functions  produce  momentum  and  energy  trans¬ 
fer  rates  that  are  zero,  in  agreement  with  the  fluctuation- 
dissipation  theorem.  The  solution  of  the  fluid  equations 
for  thermal  equilibrium  yields  the  expected  result.  See 
(65)-(69). 

The  method  we  present  has  a  formal  similarity  to  the 
standard  theory  of  anomalous  resistivity  and  anomalous 
heating. The  major  difference  is  that  in  the  work  presented 
here,  we  assume  that  the  correlation  functions  are  given  by  a 
velocity-space  Fokker-Planck  operator.  We  are  currently 
comparing  the  two  approaches  to  the  problem. 

For  nonturbulent  plasmas,  e— and  explicit  expres¬ 
sions  for  the  wave-particle  transfer  rates  were  given  for  qui¬ 
escent,  drifting,  bi-Maxwellian,  electron-ion  plasma.  Since 
the  magnetized  Balescu-Lenard  collision  operators  for  a  mul¬ 
ticonstituent  plasma  can  be  recast  into  Fokker-Planck  forms 
(see  Ref.  12),  the  effect  of  charged-particle  collisions  is  con¬ 
tained  in  the  expressions  for  and  given  by 

(52)-(58)  for  the  quiescent  plasma.  However,  if  the  plasma  is 
turbulent,  then  we  need  either  a  theoretical  or  an  experimen¬ 
tal  determination  of  e  and  {\SE^\).  As  we  mentioned  in  the 
Introduction,  if  a  solution  for  the  renormalized  propagator 
for  the  turbulent  plasma  could  be  found,  then  it  could  be 
used  to  find  the  renormalized  dielectric  screening  function 
and  the  renormalized  spectral  density  for  the  longitudinal 
electric-field  fluctuations.  These  formulas  could  then  be  used 
in  (48)-(50). 
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APPENDIX  A:  THE  VLASOV-MAXWELL  EQUATIONS 

In  this  appendix,  we  assemble  some  well-known  results 
from  plasma  theory  for  the  reader’s  convenience.  The  mate¬ 
rial  is  discussed  in  Tidman  and  Krall,^"*  Chap.  1 ,  and  we  use 
their  notation  in  this  paper. 

We  wish  to  average  the  Vlasov-Maxwell  equations  for  a 
plasma  by  splitting  the  one-particle  distribution  functions 
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and  field  quantities  into  an  ensemble-averaged  part  and  a 
fluctuating  part:  /a=(/a)  +  ^/a»  E=(E)-i-^, 

Here,  the  ensemble  average  of  a  quantity  such  as  is  de¬ 
noted  by  (fa)  and  its  fluctuation  about  the  average  value  by 
Sfa-  Averaging  the  Vlasov-Maxwell  equations  in  the  stan¬ 
dard  way,  we  obtain 


d  d  Qa 

—  -I-  V  .  —  -I- - 

dt  dr  m„ 


Qa  i 
rria 


(E)-h-vX(B) 

c 


^  [(/■«> 


dr 

dv 


av 


--( 

cdt 


477 

- j 

c 


(Al) 

(A2) 

•  <5' ,  , 

-7<B>  =  0- 

(A3) 

c  at 

(A4) 

The  charge  density  and  current  density  are  defined  as 

I  d^vifp),  j  =  2  9/3  I  (A5) 

0  J  0  J 

Here,  (/a)=(/'a(r,v,r))  denotes  the  one-particle  distribution 
function  for  the  species  a  in  the  standard  six-dimensional 
phase  space  F,  where  r=(r,v)  and  where  (fJdV  has  the 
standard  interpretation.  C«  are  the  wave-particle  correlation 
functions  for  the  problem.  It  is  important  to  note  that,  for 
quiescent  plasma,  the  contain  the  effect  of  charged- 
particle  collisions,  and  for  turbulent  plasma,  the  effect  of 
turbulent  fluctuations.  Gaussian  units  are  used,  and  (fa)  are 
normalized  so  that  fd^v{fa)=na,  where  is  the  number 
density.  The  other  quantities  have  their  usual  meaning.  The 
fluctuations  and  the  correlation  functions  determined  from 
them  satisfy  the  familiar  Vlasov-Maxwell  hierarchy  of  ki¬ 
netic  equations,  which  we  do  not  write  here. 

The  conservation  relations  for  may  also  be  deter¬ 
mined.  They  are 

dhCo,  =  Q,  (A6) 


r  <9  1  d 

X  d^vm„\C^  +  — - — (SE  X  SB)-  — 

„  J  dt4iTc  Sr 

1  d  \ 

•  (<5E^E  -f-  <5B<5B)  -i-  •  (<5E^  -i-  <5B^)I  =  0, 

477  dr  877 


(A7) 


X  d^v-mjp'Ca+——{3E}+5B?')+  — 

a  J  I  (9f077  (9r 

•— {®X^)  =  0.  (A8) 

477 

Here,  we  have  used  standard  vector  and  tensor  dyadic  nota¬ 
tion,  and  I  is  the  unit  tensor  dyadic.  In  obtaining  these  re¬ 


sults,  Maxwell’s  equations,  standard  vector  and  tensor  iden¬ 
tities,  integration  by  parts,  and  the  asymptotic  properties  of 
(/*«)  for  large  v  have  been  used.  See  Ref.  24. 

In  order  to  obtain  the  conservation  laws  for  the  number 
density,  the  momentum  density,  and  the  total  energy  density, 
we  multiply  (Al)  by  1,  and  (1/2)wq,i;^,  respectively, 
and  integrate  over  velocity  space.  We  then  substitute  for  the 
velocity  moments  of  using  (A6)-(A8).  This  gives  the 
three  conservation  laws  for  the  low-order  velocity  moments 
of  (fa)y  (E),  and  (B)  in  terms  of  the  fluctuating  electric  and 
magnetic  fields.  We  obtain 

(9  a 

-«a+  —  •««»«=  0.  (A9) 

at  ar 


— S  2  X  <5B) 


St 


477c 


+  ^  •  2  f  d^vm„yy{fj  +  -^{SE^  +  SB^)l 
ar  a  J  O'”" 

-  —(SESE  +  SBSB)  \  =  p<E)  +  -  j  X  (B),  (AlO) 
477  c 


2  +■*’>} 

+  ^  •  1 2  J  d^vjm^v^vifJ  +  X  ®)| 

=j-(E),  (All) 

where  is  the  drift  velocity  defined  as 


APPENDIX  B:  CALCULATION  OF  M^,,  W^,,  AND 
IN  THE  GUIDING-CENTER  AND  GYROTROPIC 
APPROXIMATION 

In  this  appendix,  we  give  the  details  of  the  calculation 
leading  to  Eqs.  (48)-(50).  They  may  be  obtained  as  generali¬ 
zations  of  the  work  of  Hubbard,  Ichimaru  and 
Rosenbluth,**  and  Matsuda.*^^  Uncertainties  in  the  calculation 
by  I  and  R  have  been  discussed  and  clarified  in  Refs.  43  and 
42  and  we  follow  the  procedure  of  Ref.  42  in  this  appendix. 

From  I  and  R,  the  equation  of  motion  for  a  particle  of 
charge  q  and  mass  m  is 

4v(r)  =  -®[r(f),i]  +  — v(r)XB,  (Bl) 

dt  m  me 

where  r(t)  and  v(/)  give  the  trajectory  of  the  test  particle  in 
phase  space,  B  is  the  magnetic  field,  and  SE  is  the  fluctuating 
electric  field.  We  wish  to  generalize  the  calculation  of  I  and 
R  by  considering  B  to  be  a  weak  function  of  position.  The 
meaning  of  the  weak  variation  of  B  is  given  in  Sec.  III.  In  the 
guiding-center  approximation,  we  consider  a  velocity-space 
coordinate  system  which  slowly  changes  its  orientation  as  s 
varies  along  the  B-field  line  so  that  the  axis  is  always 
oriented  parallel  or  antiparallel  to  B.  Since  B  is  changing 
slowly  with  5,  it  may  be  considered  as  locally  uniform.  In 
this  appendix,  all  subscripts  are  suppressed  in  order  to  sim- 
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plify  the  notation.  The  trajectories  of  the  particles  are  found 
by  integrating  (Bl).  They  are 

r(0  =  r(0)  +  ^HW-v(0) 

+  ^t  dt'Hit-t')-SE[r(t'),t’l  (B2) 

o  Jq 

v(/)  =  B(r)  •  v(0)  +  -  f  dt'B{t-t')-SE[rit'),t'l  (B3) 
f»Jo 

where  fl=gB/mc  is  the  cyclotron  frequency  including  the 
sign  of  q.  Here,  H  and  B  are  dimensionless  tensors  which 
characterize  the  helical  motion  of  a  charged  particle  in  a 
magnetic  field  and  are  expressed  in  Cartesian  coordinates 
(.x,y,z)  as 


H(r)  = 


BW  = 


+  sin  Clt 

-  ( 1  -  cos  Hr) 

0 


1  -  cos  Cit  0 
sin  fir  0 

0  Hr 

+  cos  Hr  sin  fir  0 
-  sin  Cit  cos  Cit  0 
0  0  1 


(B4) 


(B5) 


The  tensor  B(r)  should  not  be  confused  with  the  vector  B, 
which  denotes  the  magnetic  field.  Following  I  and  R,  the 
increments  Ar(T)  and  Av(r)  of  the  particle  coordinates  dur¬ 
ing  the  time  interval  rare  defined  as  Ar(T)=r(r)-r(0)  and 
Av(r)  =  v(T)-v(0).  They  may  be  expressed  in  terms  of  the 
fluctuating  electric  field  by  making  a  Taylor  expansion  of 
Ar(r)  and  Av(r)  in  powers  of  the  fluctuating  electric  field, 
and  then  truncating  the  expansion  to  the  second  order.  These 
formulas  are  given  in  I  and  R.  For  Av(t),  they  obtain 


Av(r)  =  [B(r)  -  I]  •  V -I- —  (  t/r'B(r')  •  <5E[ro(r- T'),r- /] -H |  ^r'B(r') 
wJo  mBjQ 

■  di"Yi{f)  ■  <5E[ro(r-  /  -  /),r-  /  -  /]j  •  ^<5E[ro(r-  /ij 


(B6) 


r 


where  ro(r)  =  r(0)  +  (l /n)H(r)-v.  In  making  this  approxima¬ 
tion,  we  have  assumed  that  the  turbulence  is  sufficiently 
weak  so  that  the  fluctuating  electric  field  has  a  small  effect 
on  the  unperturbed  particle  orbits  denoted  by  ro(/). 

The  diffusion  tensor  in  velocity  space  is  defined  as  an 
ensemble  average  of  the  dyadic  Av(r)Av(T),  as 

(Av(r)Av(r)}  ^  . 

D(v)  =  -^ — — — (B7) 
r 

The  time  r  must  be  chosen  so  as  to  satisfy  the  condition 

T>  72,  (B8) 

where  Ti  is  defined  in  Sec.  VI  and  T2  is  discussed  below. 
Since  Hubbard*^  has  shown  that  the  self-field  contribution  to 
D  is  small  compared  to  the  dominant,  second-order-in-® 
contribution,  and  since  (®)=0,  we  may  expand  AvAv  using 
(B6),  take  the  ensemble  average,  and  keep  only  the  second- 
order  terms.  The  dominant  second-order  contribution  is 

D(v)  =  -^f  t//  f  dmT') 

•  ®[ro(T-  r'),r-  +f) 

•  <5E[ro(T- t' -r),T- / -r]>,  (B9) 

where  the  condition  (B8)  has  been  used  to  extend  the  range 
of  the  t  integration  from  -oo  to  This  is  (13)  of  I  and  R. 
The  other  two  second-order  terms  are  given  in  Ref.  42.  Us¬ 
ing  the  procedure  to  evaluate  the  t'  integration  that  we  de¬ 


scribe  below,  which  is  used  throughout  this  appendix,  it  is 
possible  to  show  that  the  other  two  second-order  terms  van¬ 
ish.  Therefore,  we  are  left  with  (B9)  for  D.  The  procedure  we 
use  to  evaluate  the  t'  integration  is  the  same  as  that  used  in 
Ref.  42  and  is  described  as  follows.  When  no  magnetic  field 
is  present,  Hubbard  noted  that  correlations  of  the  fluctuating 
electric  field  persist  only  for  times  on  the  order  of  / 

When  a  magnetic  field  is  present,  Matsuda"^^  has  assumed 
that  the  major  contribution  to  the  t'  integration  comes  from 
times  where  0^r'~  a  few  times  72,  where  72(=7^(.)  is 
the  autocorrelation  time  for  the  turbulence  and  where 
We  accept  this  assumption  and  will  evaluate  be¬ 
low  the  t'  integration  in  (B9),  and  the  analogous  integration 
in  also  given  below,  by  integrating  on  t'  and  then  taking 
the  limit  as 

We  now  generalize  the  I  and  R  calculation  to  include  a 
weakly  varying,  fluctuating  electric  field  in  space  and  a 
slowly  varying,  fluctuating  electric  field  in  time  by  introduc¬ 
ing  the  tensorial  spectral  density  of  the  electric-field  fluctua¬ 
tions 

<®(r,r)®(r',r')) 


Xexp[/k  •  (r-  r')  -  (BIO) 

Here  (®(r,r)®(r' ,r'))  is  the  two-point,  two-time  correla- 
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tion  function  for  the  turbulent  electric  field.  It  depends  not 
only  on  r-r'  and  t-t\  but  on  the  average  spatial  (r-i-r')/2, 
and  average  temporal  coordinates.  In  order  to  re¬ 

mind  us  that  we  are  describing  correlation  functions  for  the 
turbulent  plasma,  we  use  the  tilde  notation.  Now,  we  wish  to 
impose  the  same  restrictions  on  the  length  and  time  scales  of 
the  average  spatial  and  the  average  temporal  variations  of 
{^®*[(r-i-r')/2,(r-i-r')/2;k,cu])  as  we  did  in  Sec.  Ill  for 
the  variation  of  (E).  Therefore,  we  assume  that  the  length 
scale  for  the  variation  of  with  (r-i-r')/2  in  the  di¬ 

rection  perpendicular  and  parallel  to  B  is  large  compared  to 

and  an,  respectively.  We  also  assume  that  the  time  scale 
for  the  variation  of  {^^*)  with  (r-i-r')/2  is  large  compared 
to  |27r/n|.  In  the  guiding-center  coordinate  system  we  use, 
which  was  discussed  in  Sec.  Ill,  {r+v')l2-^s  and  {t-\-t')l2 
— >r.  Here,  s  and  t  are  the  weak  spatial  and  slow  temporal 
coordinates,  and  r-r'  and  t-t'  are  the  strong  spatial  and  fast 
temporal  coordinates,  where  5  is  the  distance  along  the  mag¬ 
netic  field  line.  Again,  it  is  important  to  point  out  that  the 
scaling  assumptions  leading  to  (BIO)  for  the  spectral  density 
are  valid  only  for  a  large  separation  of  the  length  and  time 
scales. 

When  the  fluctuations  are  electrostatic,  the  tensor  of  the 
spectral  density  of  the  electric-field  fluctuations  has  a  simpli¬ 
fied  form, 

(®^*(j,r;k,cu))  =  (kk/A:")<|^P|(5,r;k,cu)>,  (B1 1) 


Xexp{/(A2n-i-^||i;||-  cu)r}.  (B15) 

In  writing  (B15),  we  have  chosen  a  Cartesian,  velocity-space 
coordinate  system  where  e^,  e^,  and  are  appropriate  unit 
vectors  and  the  magnetic  field  is  oriented  either  parallel  or 
antiparallel  to  e,.  Here,  <p(^)  is  the  angle  that  the  component 
of  v(k)  in  the  e^-e^  plane  makes  with  the  axis.  Also, 
where  is  the  component  of  v(k)  in  the 

Cj^-Cv  plane,  and  is  the  component  of  v(k)  either  par¬ 
allel  or  antiparallel  to  the  magnetic  field.  Using  (B5)  and 
integration  by  parts,  we  see  that  the  r'  integration  can  be 
carried  out  and  the  limit  as  nr/27r— ^0  taken,  as  described 
above,  to  give 

2  r  ci^lc 
m^J  (27r)^J_oo 

r"  1 

X  ^y/-[kB(r)  •  k]exp[-  ik  •  r)  •  v  -  icot]. 

J-oo  k- 

(B16) 

This  is  Eq.  (17)  of  Ref.  42  when  is  independent  of  s 

and  t.  It  is  the  generalized  D  when  has  a  weak  spatial 

and  a  slow  temporal  dependence  before  the  gyrophase  aver¬ 
age  is  calculated.  From  Sec.  VI,  we  see  that  since  /„  is  re¬ 
placed  by  7„,  we  need  only  to  calculate  the  gyrophase  aver¬ 
age  of  D  in  order  to  evaluate  and  Therefore, 


where  (|^£'^|(5,r;k,co))  is  the  spectral  density  of  the  turbu¬ 
lent,  longitudinal  electric-field  fluctuations  and  has  the  sym¬ 
metry  (|<SE^|(5,/;k, (o))=(|^E^|(5,r;-k,-a>)).  Using  (BIO), 
we  find  that  the  generalized  diffusion  tensor  is 
2  r  i3i  r+oo 


X 


iyi/4i 


[B(r')-k][B(T'+f)-k] 


Xexp[/k  •  Aro(/)  -  iwf], 


(B12) 


Aro(f)  =  ro(r-  /)  -  ro(r-  r*  -  0 

=  a-'[H(T-T')-H(T-7'-r)]-v.  (B13) 

Here,  the  slow  time  dependence  of  {|<5£^|)  has  been  excluded 
from  the  integration  over  the  fast  time  scale.  Using  the 
identity^* 

+0C 

exp(/z  sin  i/f)  =  2  (B14) 


-  f  (p’k  /  I 


X 


k\  cos  Or  -  k\  sin  Clt  0 
k\  sin  fit  +  k\  cos  Clt  0 


0  2k\ 

X  2  JlU)t^p[i{nfl  -I-  k\\Vii  -  (o)t] 


0 

+00 


(B17) 


Noting  that  only  the  even  part  on  /  survives  the  inversion  on 
k,  (i),  and  n,  we  obtain 


D  = 


0 

0 


0 

5 

0 


0 

0 

A, 


where 


(B18) 


we  find  that  the  expanded  expression  for  D  is 

2  f  j3i.  "=+“'=+* 

I  77”^  2  2  ^fw<|<5p(^,^;k,a))|> 

W  J  (277)^ 

X  j  d7'-^[B(7')  ■  k][B(T'  +  /)  •  k] 


X(|<5E"|(5,r;k,a>))^nri^-i-  -  co), 

7,2 


(B19) 


27rq^  f  (pk  f'*’*  V'  r  2  /  x  ^  /  xn 

X(|5£2l(5,r;k,  w))^nn  +  -  w).  (B20) 


X  y„(fV/(^)exp{j(n  -  l)[<p-  dj- CI(t-  t')]} 


Equations  (B19)  and  (B20)  are  the  generalizations  of  Eqs. 
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(26)  and  (25)  of  I  and  R  to  include  the  weak  spatial  and  slow 
temporal  variation  of 

In  order  to  find  the  generalized  friction  coefficient,  again 
we  follow  I  and  R.  It  is  defined  as 


F(v)  = 


(Av(t)) 

r 


(B21) 


where  r  is  a  time  interval  satisfying  (B8).  As  shown  by 
Hubbard,*^  there  are  two  contributions  to  F:  a  part  due  to  the 
correlated  effect  of  the  electric-field  fluctuations  on  the  par¬ 
ticle  orbit,  F^,  and  a  polarization  part  due  to  the  self-field  of 
a  test  particle  in  a  plasma,  F^, 

¥  =  F^+¥P.  (B22) 

The  formulas  given  in  I  and  R  may  be  generalized  in  the 
same  way  as  we  did  to  obtain  (B 1 6)  for  a  weakly  varying  B 
field  with  position,  and  a  weakly  and  slowly  varying  {\SE'^\) 
with  position  and  time,  respectively.  In  the  calculation  of  F^, 
we  note  that  the  first  two  terms  in  (B6)  do  not  contribute 
because  their  ensemble  average  is  zero,  and  that  the  contri¬ 
bution  comes  from  the  third  term.  Following  procedures 
similar  to  the  calculation  of  D,  we  obtain 

Xexp[-  ik  •  n~'H(-  t)  icot].  (B23) 

We  have  used  the  above-described  procedure  to  evaluate  the 
t'  integration  and  (B8)  to  extend  the  t  integration,  first  from 
0  to  00,  and  then  from  -oo  to  +oo  by  taking  half  the  value. 
Equation  (B23)  is  the  same  as  Eq.  (24)  of  Ref.  42  when 
(|^£^|)  is  independent  of  s  and  t  and  is  the  generalized  F^ 
before  the  gyrophase  average  of  v  F^  is  calculated.  It  is  also 
possible  to  show  that 

d/^’D  =  2F^  (B24) 


where  D  is  given  by  (B16). 

The  second  term  in  (B22)  is  due  to  the  polarization  of 
the  plasma  and  has  been  given  when  there  is  no  magnetic 
field  present  by  Hubbard.*^  When  a  magnetic  field  is  present, 
we  obtain 

[  77^3  f  &>)■'] 

m  J  (277)^  J_„ 

xj  <fr|^jexp[-;k-n-'H(/)-v  +  /a)/].  (B25) 

Here,  e  is  the  dielectric  screening  function  for  the  turbulent 
plasma.  It  varies  weakly  with  s  and  slowly  with  t  as  the 
properties  of  the  plasma  medium  vary.  We  also  note  that 
E(5,/;k,a>)*  =  £(5,r;-k,-a>),  and  (B25)  is  the  same  as  (30) 
of  I  and  R  and  (27)  of  Ref.  42  when  e  is  independent  of  s 
and  t. 


From  Sec.  VI,  we  see  that  in  order  to  find  VVn  and 
we  need  to  find  the  gyrophase  average  of 


V  •  (F^+  F'’)/=  {v  •  (F/+  F'’)},/+  {v  •  (F/+  F")} 


(B26) 


We  proceed  as  in  the  calculation  for  D  to  obtain 

rlir 


r 

(27r)  ^  I  d(pm(\  •  F^)| 

Jo 

_q^C  r 
m  J  (277)^ 


X  2  [y„'.,(a-4i(^)]<l<5£'l(^./;k,a,)) 


X  S{nCl  -H  A:||i;||  -  w). 


(B27) 


(277)  ^  I  d(pm(\  •  F^)j 


■'/ 

Jo 

_  ^  f  n 

m  J  (277)0- 


X  2  [y„^,(a+-/„^i(l)](|<5p|(5,nk,a>))  + 
/!=— 00  \ 


n=-x 

+x 


<7tl> 


(B28) 


A  similar  calculation  for  F'’  yields 

r2ir 


r2TT 

(277)-' J  d^miy-YO), 


X  2  y^(^)^nn  +  A:||ii||- w), 


d(pm{\  •  F^)  J 


] 

(B29) 


70 

g^TT  f  d^k  r  /  N-1 


1 


x2  nflJ^(^)S(nfl  +  kiiVii- oj).  (B30) 

n=-x 

Combining  (B27)-(B30)  with  (B19)  and  (B20)  allows  us  to 
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write  expressions  for  and  in  the  guiding-center  and 
gyrotropic  approximation.  They  are  given  by  (49)  and  (50). 
We  note  that  and  the  first  term  on  the  right-hand  side  of 
(B28)  cancel,  and  that  does  not  appear  in  (50). 

In  order  to  find  we  need  to  evaluate  the  gyrophase 
average  of  and  Using  similar  procedures,  we  obtain 


X(|^£^|(5,/;k,co)}  r  S(nCl  -i-  k^Vw  -  cu),  (B31) 


Xlm[e(5,/;k,ct;)  -l- /:||i;||  -  cu).  (B32) 


Combining  the  two,  we  obtain  (48). 
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